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755. 


ON THE MATRIX (a, 0) AND IN CONNEXION THEREWITH 
OK +8 


THE FUNCTION +? 


cx+ d 


[From the Messenger of Mathematics, vol. 1x. (1880), pp. 104—109.] 


In the preceding paper, [due to Prof. W. W. Johnson,] the theory of the symbolic 


powers and roots of the function cote is developed in a complete and satisfactory 


manner; the results in the main agreeing with those obtained in the original memoir, 
Babbage, “On Trigonometrical Series,” Memoirs of the Analytical Society (1813), Note 1. 
pp. 47—50, and which are to some extent reproduced in my “Memoir on the Theory 
of Matrices,” Phil. Trans., t. CXLVIII. (1858), pp. 17—37, [152]. I had recently 
occasion to reconsider the question, and have obtained for the nth function $x, where 
px =n , a form which, although substantially identical with Babbage’s, is a more 
compact and convenient one; viz. taking A to be determined by the quadric equation 


(A+1)?_ (a+d) 
»~  £ad—be’ 


the form is 
gn (VaR) Getta eS b) 
~ (A"*1—1) (cv +d)+(A"—2)( cs — a)’ 
The question is, in effect, that of the determination of the nth power of the 
matrix ( a, b ); viz. in the notation of matrices 
| c, d 
(2, y) = ( a, b ) (æ, y), 
c, d | 
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means the two equations m, =a%s + by, y,=cx+dy; and then if a, y, are derived in 
like manner from æ, yı, that is, if 2=a2,+by,, Y= ce, +dy and so on, En, Yn will 
be linear functions of æ, y; say we have a,=Q,0+ ny, Yn=Cnt+dny: and the nth 
power of ( a, b ) is, in fact, the matrix ( dy, bn ). 


c, d | | Cn, Un | 
In particular, we have 
(a,b F, =( @, b ), =( a®+be, b(a+d) ), 
E PASI | c(a+d), d+be | 
and hence the identity 
(a, b P—(a+d)( a, b )+(ad— me Tee yea Os 


| «, d | G: d | |0, 1 | 
viz. this means that the matrix 
( a—(a+d)at+ad—be, b,-—(a+d)b j= ('O; 0"); 
Mn ere . da(a+dd+ad—te| |0, 0| 


or, what is the same thing, that each term of the left-hand matrix is =0; which is 
at once verified by substituting for æ», bz, Cə, dą their foregoing values. 


The explanation just given will make the notation intelligible and show in a 
general way how a matrix may be worked in like manner with a single quantity: 
the theory is more fully developed in my Memoir above referred to. I proceed 
with the solution in the algorithm of matrices. Writing for ‘shortness M= ( a; b ), 

| c, a | 
the identity is 
M -— (a+d) M + (ad — bc) =0, 


the matrix ( 1, 0O ) being in the theory regarded as =1; viz. M is determined by 
CAA 
a quadric equation; and we have consequently M”=a linear function of M. Writing 
this in the form 
M*—-AM+B=0, 
the unknown coefficients A, B can be at once obtained in terms of a, B, the roots 


of the equation 
w—(a+d)u+ad—be=0, 


viz. we have 
—Aa+B=0, 


-AR +B=0; 


or more simply from these equations, and the equation for M”, eliminating a, 8, we 


have ` : 
M”, M, 1|=0; 


5 egy 


B", 8, 1 
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that is, 
M” (a— 8) —M(a" — B”) + aß (a= — Bn) = 0. 


But instead of a, 8, it is convenient to introduce the ratio à of the two roots, 


say we have a=A8; we thence find 


(A+1)8=a+d, 


rAB? = ad — be, 
giving 
(A+ 1? (a+dp 
- yr ad—be 
for the determination of A, and then 
a+d 
B br x Ae 1 ? 
het (a+d)x 
> El 


The equation thus becomes 


M"(x—1) B— M(A*®—-1) B+ A-E =0, 
or we have 
un = o- 1) M —(A"—2) B}. 


It is convenient to multiply the numerator and denominator by A+ 1, viz. 


thus have 


m= Ê [anh —1) M+) [(M-O+1) 8)) 


The exterior factor is here 
ee at sy 
A SI E TEE | : 
moreover (A+1)8 is =a+d: hence 
Ma C ae 6), 
le, al 


and 
M-—(.+1)B8=( a, b )—(at+d, 0 D), =(-d b_); 
mie | Oo awa aed 
the formula thus is 


aes sai es G P —1) i ‘ f y+ A-A) i ie a (i ; 


viz. we have thus the values of the several terms of the nth matrix 
M n an» Da Jz 


Cn, dy, 
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and, if instead of these we consider the combinations a,#+b, and c,v+d,, we then 
obtain 


SA) (Om = 1) (a+b) +O") (detb), 


Crt+dn= , (a — 1) (ca+d)+(A"—A)( cx—a)}; 


1 
dnt + by = = — ( 


and in dividing the first of these by the second, the exterior factor disappears. 

It is to be remarked that, if n=0, the formule become as they should do œs +b, =z, 
cs +d=1; and if n=1, they become aav+b,=ax+b, x+d,=cx+d. 

If »¥”™—1=0, where m, the least exponent for which this equation is satisfied, is 


for the moment taken to be greater than 2, the terms in { } are 


i (\—1)(ax+b)+(1—)(— de +b), 
an 
A—1)(cæ +d)+(1=A)( ce —a); 


viz. these are (à—1)(a+d)x, and (A—1)(a+d), or if for (A—1)(a+d) we write 


a+d 
(2 a erni , the formule become for n =m 
a + d\” 
AmE + bm = F) a, 
Cmt +d Ten 
MY r aa oa TE 1 , 
viz. we have here 
AmE + bm 
Gra + dn 


or the function is periodic of the mth order. Writing for shortness ee s being 
n 


any integer not =0, and prime to n, we have A= cos 29 +7sin 29, hence 
1 += 2 cos ÙX (cos $+7sin X), 


2 . 
r A =4 cos; consequently, in order to the function being periodic of the nth 
order, the relation between the coefficients is 
sir _(@ +d 
iii n ad—be’ 


The formula extends to the case m= 2, viz. cos 4(sm)=0, or the condition is 
a+d=0. But here ~-+-1=0, and the case requires to be separately verified. Recurring 
to the original expression for M?, we see that, for a+d=0, this becomes 


a+b, 0 |, =(a?+6c)| 1, 0 ? 
i 0 , d&+be (0 ane I 
that is, 
dye + by _ 
u Fila F 


or the result is thus verified. 
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But the case m=1 is a very remarkable one; we have here X=1, and the 
relation between the coefficients is thus (a+d)*=4(ad—bc), or what is the same thing 
(a—dy+4be=0. And then determining the values for X=1 of the vanishing fractions 
which enter into the formule, we find 


1 


ga (a + dy {(n + 1) (ax +b) + (n — 1) (— dæ + by}, 


Ant + by = 


ose + dy = =; (a + d) (x +1) (ca +d) +(n—1)( cv — a)}, 
or as these may also be written 
Ant + bn = = (a +d)" {a [n (a — d) + (a + d)] + 2nb}, 


Bs 


Cnt + dn, i Qn 


(a +d)" {æ . 2nc +[—n(a-—d)+a+d}]}, 

which for »=0, become as they should do a+ b =<, qa+d,=1, and for n=1 they 
a,x + b, 
ce + d 
function is not periodic of any order. This remarkable case is noticed by Mr Moulton 
in his edition (2nd edition, 1872) of Boole’s Finite Differences. 


become a,¢+b,=axr+b, oa+d,=ce+d. We thus do not have =a, and the 


If to satisfy the given relation (a — d} + 4bc = 0, we write 2b = k (a — d), 2c = — : (a—d). 


then the function of a is 
ax+4k(a—d) 
—tk>(a-—d)a+d’ 


and the formule for the nth function are 


Git + bn = z; (0 +A) (a+) e +n (a-d) (+h), 


ent +da = 5, (a + d) fara —n(a-d)(F+ 1)}; 
which may be verified successively for the different values of n. 


Reverting to the general case, suppose n=, and let w be the value of $% (a). 
Supposing that the modulus of à is not =1, we have A” indefinitely large or 
indefinitely small. In the former case, we obtain 


_dA(awt+b)+(—da+b) —_(Aa—d)e#+bA+l1), 
"=X (ew +d)+( ce—a)’  c(rA+1)a+Ad-a’ 


which, observing that the equation in A may be written 


Aa—d _b(A+1) 
c(A+1) Ad-a’ 
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is independent of æ, and equal to either of these equal quantities; and if from these 
two values of u we eliminate A, we obtain for u the quadric equation 


cu? —(a—d)u—b=0, 
that is, 
_autb 
KAE 


as is, in fact, obvious from the consideration that n being indefinitely large the nth and 
(n+1)th functions must be equal to each other. In the latter case, as A” is indefinitely 
small, we have the like formule, and we obtain for w the same quadric equation: 
the two values of u are however not the same, but (as is easily shown) their product 
is =—b+c; wu is therefore the other root of the quadric equation. Hence, as n 
increases, the function $2 continually approximates to one or the other of the roots 
of this quadric equation. The equation has equal roots if (a—d)?+4bc=0, which is 
— 2b 
a—d’ 
which result is also given by the formule of the special case on writing therein n= œ. 


s A ‘ } : 1 
the relation existing in the above-mentioned special case; and here w= 5 (a—d), = 


C. XL 33 
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